Triangulating tunneling resonances in a point contact 
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We observe resonant tunneling in silicon split gate point contacts implanted with antimony and 
defined in a self-aligned poly-silicon double gate enhancement mode Si-MOS device structure. We 
identify which resonances are likely candidates for transport through the antimony donor as op- 
posed to unintentional disorder induced potentials using capacitance triangulation. We determine 
the capacitances from the resonant feature to each of the conducting gates and the source/drain two 
dimensional electron gas regions. In our device geometry, these capacitances provide information 
about the resonance location in three dimensions. Semi-classical electrostatic simulations of capac- 
itance, already used to map quantum dot size and position, identify a combination of location and 
confinement potential size that satisfy our experimental observations. The sensitivity of simulation 
to position and size allow us to triangulate possible locations of the resonant level with nanometer 
resolution. We discuss our results and how they may apply to resonant tunneling through a single 
donor. 
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I. INTRODUCTION 

Numerous groups are developing single electron spin 
devices in silicoiP^. Electron spins in silicon can have 
very long coherence times^ ^, which is desirable for both 
the basic study of spin interactions using engineered few 
spin systems as well as the pursuit of applications such as 
quantum computing and spintronics. Donor atoms pro- 
vide a natural confinement potential for single electrons 
at low temperature. Donors placed intentionally (i.e., im- 
plant) or adventitiously into tunnel barriers of silicon de- 
vices have been probed using transport spectroscopy and 
charge sensing^ ^ ^ . Single electron spin measurements 
have been reported recently in intentionally implanted 
structures introducing a tantalizing fabrication path to- 
wards future single donor spin devices^ ^. Two common 
challenges to developing single donor devices are estab- 
lishing strict control over the position of the donor and 
avoiding interference from defects. 

In this letter we report fabrication, measurement and 
modeling of a self-aligned, donor-implanted nanostruc- 
ture that provides a path towards an improved alignment 
of donors to neighboring ohmic leads and lateral charge 
sensor structures. Measurements of implanted devices 
show more tunneling resonances in implanted tunnel bar- 
riers, with respect to unimplanted tunnel barriers. The 
resonant tunneling locations are triangulated using mea- 
sured capacitance ratios. Triangulation reveals that some 
resonances likely correspond to transport through regions 
that are not intentionally implanted. This highlights 
the value of geometrical analysis, such as triangulation, 
to improve the confidence in identification of resonances 
that are candidates for transport through intentionally 
implanted donors instead of defects. The technique of 
triangulation is, furthermore, more generally applicable 
to any system for which the presence of defects intro- 



duces ambiguity about the source/location of single elec- 
tron behavior, for example, silicon, graphene and carbon 
nanotubes. 



II. EXPERIMENT 

We fabricated and measured a point contact formed 
by a split gate geometry in an enhancement mode Si- 
MOS nanostructure. Fig. [ija). Poly-silicon (poly), de- 
posited 200 nm thick and then degenerately doped via ion 
implantation, form first level gates, used to deplete the 
electrons and form the point contact. They are isolated 
from the silicon substrate by a 35 nm thermal Si02 layer. 
The second level, an enhancement gate, is formed with a 
100 nm aluminum layer, covering the entire active region 
of the device and overlapping highly doped substrate re- 
gions which serve as ohmic contacts to the induced elec- 
tron layer. The Al top gate is isolated from the poly by 
a 60 nm aluminum oxide layer grown by atomic layer de- 
position (ALD). Quantum dots with periodic Coulomb 
blockade and lateral charge sensing of the quantum dot 
electron occupation have been demonstrated with similar 
processing^ ^. 

A small number of antimony atoms were implanted in 
to one of two point contacts in the device. After the poly 
gate pattern is etched. Fig. [ijb), the device was covered 
with 300 nm of PMMA, and an 80x80 nm^ opening was 
formed in one constriction (left side, blue arrow), while 
a second constriction remained unexposed (right side, 
brown arrow). A blanket 100 keV Sb implant with dose 
of 4 X 10^^/ cm^ was used. SRIM simulations^^ project 
that the combined thickness of the 35 nm gate oxide and 
the PMMA should result in Sb concentrations below the 
p-type background of the substrate, > 100 Ohm-cm or 
~ 10^^ /cm^. The dose and energy correspond to an av- 



erage of 8 Sb ions in the 80x80 nm^ window, at a depth 
of 10 nm below the Si/Si02 interface, according to sim- 
ulations done with SRIM. Roughly 50% of the dose is 
deposited in the Si02 dielectric. The implanted constric- 
tion acts as a control case, to characterize non-implant 
process damage that might mimic donor assisted tunnel- 
ing through the split gate barrier. 

The device structure is completed by reoxidation at 
900 °C for 24 minutes followed by the addition of ALD 
AI2O3 and Al metallization, as described previously^. 
The reoxidation step activates the dopant and anneals 
damage. High frequency and quasi-static capacitance- 
voltage measurements were done to characterize the 
residual damage for a similar process flow using blan- 
ket depositions, implants and large Al capacitor struc- 
tures, shown in Table [ij To within the uncertainty of 
the measurement, reoxidation removes residual damage 
from implant whereas a rapid thermal anneal leaves non- 
negligible defect charge in the oxide, which could intro- 
duce a higher number of defect related resonances. A 
simple estimate for the diffusion length of Sb for this 
thermal budget is 6 nm^^. We note that because Sb is 
a vacancy based diffuser, excess interstitials introduced 
due to implant and oxidation that enhance other dopant 
diffusion will not greatly enhance Sb diffusion, and may 
suppress itP^. 

In Fig. [ijc), we plot a schematic model of the elec- 
tron potential in the Si channel, along the blue curve in 
part (b). At the extremes are high-density source and 
drain regions, where the Fermi energy is well above the 
conduction band edge. Between the two is a depletion 
region, where the negative bias on the depletion gates 
has brought the conduction band edge above the Fermi 
energy, and brought an energy level of a hypothetical 
donor or trap below the band edge into resonance with 
the source and drain. 

We examine the point contacts using transport spec- 
troscopy. The threshold voltage through the implanted 
point contact is appreciably higher than the control, 17 
and 10 volts, respectively. To determine the threshold, 
we extrapolate to zero conductance from a region where 
the conductance is linearly related to the enhancement 
gate bias, fl] (d) shows the current through the implanted 
point contact (left side, blue arrow in Fig. [ijb)) as a 
function of depletion gate bias. We measure the conduc- 
tance using low- frequency lock-in techniques, in this case 
a 15 — 100/i Vrms AC excitation at 27 Hz, and measured 
the resulting current. 

We observed numerous peaks in the sub-threshold con- 
ductance, in contrast to a single small peak observed in 
the control side. Several other samples show a similar 
trend of increased numbers of resonances on the Sb im- 
planted point contact compared to controls. A point im- 
planted with silicon, dose 4 x 10^^/ cm^ and energy 45 
keV, chosen for its similar distribution of Si interstitials 
as the Sb implants, did not show a substantial increase 
in resonances below threshold. Sub-threshold conduc- 
tance resonances through a point contact normally result 
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FIG. 1: (color online) (a) Cross section of our device, not to 
scale, showing polysilicon depletion gates and aluminum in- 
version gate, (b) Scanning electron micrograph of the polysil- 
icon depletion gate layout, from a representative sample, (c) 
Schematic energy diagram, along the blue arrow in part (b). 
(d) Conductance through the left constriction (blue arrow in 
part (b)) at T = 20 mK, showing sub-threshold tunneling 
resonances (labeled A, B, and W). 



TABLE I: Device fabrication process damage parameters, 
from capacitance vs. voltage measurements on planar test 
devices. 



Sample type 


Qfb 
(q/cm^) 


D., 

(q/cm^eV) 


Normal process 


+5 X 10^° 


1 X 10^° 


Silicon implant + 30 sec @ 800 C 
(Donor device equivalent dose) 


+3 X 10^^ 


2 X 10^^ 


Silicon implant + 24 min @ 900 C 


+ 1 X 10^^ 


1 X 10^° 



from tunneling through a localized confinement potential 
within the tunnel barrier. The confinement can be due to 
intentionally implanted donor potentials or unintentional 
damage related defects (e.g., charged defects in the oxide) 
or impurities (e.g., metals or other background dopants). 
We will discuss results primarily from the resonances la- 
beled A, B, and W. 

We measured the tunneling current dependence on 
each gate bias, as well as the source-drain bias across the 
constriction. Several sets of data, each taken with a dif- 
ferent configuration, are shown in Fig. |2|b - d). Peaks A, 
B, and W (not shown here) all exhibit clearly defined tri- 
angular regions of high conductance consistent with dis- 
crete electron. Coulomb blockaded transport. If the drain 
2DEG region is held at ground, the slope of the positive- 
going line demarcating the high-conductance region is 
a measure of the ratio of capacitances Cgate/{C — Cs)- 
Here C is the total capacitance from the localized elec- 
tron state to all conductors in the device, Cs is the ca- 
pacitance between the source 2DEG region and the lo- 
calized state, and Cgate is the capacitance between the 
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FIG. 2: (a) Close up SEM image of the implanted constric- 
tion, with labels for the depletion gates used in this analy- 
sis, (b) Conductance as a function of depletion gate bias and 
top (aluminum) gate bias. The movement of peaks A and 
B, dVwire/dVrG, is determined by the ratio of capacitances 
Ctop/Cwire- (c) Conductance vs. source-drain bias and gate 
bias (both qpc and wire gates). Source grounded, (d) Con- 
ductance vs. source-drain bias and gate bias (wire gate only, 
qpc fixed). Drain grounded. 



depletion gate and the localized state. The slope of the 
negative-going line is a measure of the ratio Cgate/Cs- 
By performing this measurement for both the wire and 
qpc gates, and by switching which 2DEG region is held at 
ground potential, we obtained a set of capacitance ratios 
for each peak. A set of measurements for the constriction 
determines Cs/C, Cd/C, C^ire/C, Cqpc/C, and Ctop/C, 
which together account for more than 97% of the total 
capacitance (Fig. |2|c) and (d)). To determine the top 
(aluminum) gate capacitance ratio, we measure the po- 
sition of each peak as a function of top gate and (e.g.) 
wire gate voltage, which determines the ratio Ctop/Cwire- 
Note that for the measurements of peak W, an incom- 
plete set of data allows us to only determine the total 
Cwire + Cqpc/C^ but wc cauuot determine separate val- 
ues for Cwire/C or Cqpc/C. 

We note that the measured capacitance ratios are 
found to be self-consistent. That is, a capacitance ratio 
such as Cd/C can be directly measured or indirectly esti- 
mated through measurement of the ratios Cs/C^ Cqpc/C^ 
Cwire/ C^ and Ctop/C^ in conjunction with the constraint 
that the sum of all capacitance ratios must equal unity. 
All methods of determining each ratio are in good agree- 
ment with one another to within 5%, despite uncertainty 
due to drift from measurement to measurement, and dif- 
ficulty in precision determination of the slope within a 
single measurement. 

Assuming the capacitance is primarily a geometric ef- 
fect, these ratios contain information about the spatial 
location of the bound electrons. Different ratios suggest 
different spatial position or size of the confined electron 
at each resonance. Measuring Cwire /C and Cqpc/C de- 
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FIG. 3: (a) An oblique view of the 3D structure of our device, 
used for calculating the 2D electron density in the channel, 
shown in (b). (b) Constant electron density contours, from 
1 X lO^V cm^ to 1 X lO^V cm^ in steps of 1 x lO^V cm^ We 
take the size of the source and drain regions to be that area 
with density higher than 1.5 x 10^^/ cm^. 



termines the position of the localized state perpendicular 
to the current direction, while Cs/C and Cd/C deter- 
mine the position along the current direction. We have 
tabulated the results for each of the three points of in- 
terest in Table [n 



III. MODELING 

The measured capacitance ratios bound the possible lo- 
cation of a single electron on a donor through the geomet- 
ric dependences of the capacitance to position. A detailed 
capacitance model is used to establish quantitatively the 
predicted position (i.e., capacitance triangulation) and 
confidence in the triangulation. A semi-classical Thomas- 
Fermi calculation done in TCAD Sentaurus^^ is used to 
establish the position and density of the 2DEG. Then ca- 
pacitances are calculated using the commercial finite ele- 
ment modeling package CFD-ACE^^^, in which all poly 
gates, the top Al gate and the S/D 2DEGs are modeled 
as perfect conductors. A 2DEG density of 1.5 x 10^^/ 
cm^ is used as the metallic edge for the 2DEG. Any den- 
sity above 1.5 x 10^^/ cm^ we regard as perfectly con- 
ducting, and below as perfectly insulating. This value is 
similar to an experimentally measured metal-insulator- 
transition temperature for MOS 2DEG&^^ and has been 
successfully used in this way for modeling quantum dot 
capacitances^^. The capacitances between the gates and 
a 1.8 nm radius sphere, the Bohr radius of the first bound 
electron on antimony in bulk silicon, are calculated in this 



TABLE II: Measured vs. simulated capacitance ratios for three tunnel resonances in our device. 





B 
Expt. 


Sim. 


A 
Expt. 


Sim. 


W 
Expt. 


Sim. 


CsIC 


0.436 ±0.018 


0.435 


0.378 ±0.014 


0.379 


0.257 ±0.024 


0.257 


Cd/C 


0.393 ±0.010 


0.392 


0.474 ± 0.006 


0.473 


0.627 ±0.053 


0.626 


^wire/ ^ 


0.014 ±0.013 


0.013 


0.019 ±0.001 


0.019 


0.094 ±0.009 


0.0025 


Cqpc/C 


0.134 ±0.014 


0.132 


0.099 ±0.007 


0.095 


0.060 
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and drain regions, Pt. B has Cs ^% larger than C^, Pt. 
A has Cd 11% larger than Cs] the resulting perpendicu- 
lar positional shift (from source to drain) for Pt. B is ±2 
nm and —4 nm for Pt. A. This highlights the potential 
for high precision triangulation using this method, but 
also highlights the potential for error due to uncertain- 
ties in the exact gate dimensions. Lateral dimensions of 
the poly are measured by scanning electron microscopy 
on a representative sample. We estimate that the uncer- 
tainty in dimensions is dominated by sample to sample 
variation, and is less than or equal to 15 nm for gaps 
and widths. 



FIG. 4: Best fit locations for resonances A, B, and W. The 
simulation determines a unique position for each possible com- 
bination of measurements. These error bars reflect exper- 
imental error, and do not include uncertainty from possible 
disagreement between the model and the actual device, which 
is expected to be small. 



case to examine what antimony locations would be con- 
sistent with the measurements. IsFa) is a view of the 3D 
model used for this calculation, and a contour plot of elec- 
tron density for the gate biases used in the measurements 
of[2|b) is shown in|3jb). A diagram of all the conductors 
in the simulation, labeled CI to C19, is shown inlsFc). 
Conductors C3 and C4 represent a single trial position for 
a donor in the tunnel barrier regions. We chose positions 
by trial and error until a good match to the measured 
capacitance ratios was found. Further refinement of the 
position was established through interpolation between 
the database of previous positions and capacitances, |4] 
The interpolation produces an optimal single position for 
each resonance. Negative values indicate a position under 
the poly gate (i.e., an unlikely location for an intention- 
ally implanted Sb). The W peak position, triangulated 
by capacitance ratios, is the only resonance that is consis- 
tent with transport through a location between the poly 
split gates where Sb was targeted. 

The capacitance model indicates that small differences 
in capacitance ratios can correspond to significant po- 
sitional shifts of many nanometers. For example, the 
measured capacitances ratios {CwirelC and Cqpc/C) for 
Pts. A and B differ by 24-36%; the resulting lateral po- 
sitional shift (along the QPC gate) from Pt. B to Pt. A 
is 54 nm as determined by simulation, [4J For the source 



IV. DISCUSSION 

In three dimensions, we can determine lateral position 
relatively easily (position parallel to the silicon-oxide in- 
terface), but vertical displacement is less direct. The 
capacitance from source or drain to the localized state 
determines both the length of the tunneling region that 
forms in the constriction, as well as any displacement 
from the center of the tunnel barrier. If we assume that 
the source and drain regions have similar shapes, any dif- 
ference in source vs. drain capacitance must be due to 
displacement from the mid-point between the source and 
drain in the tunnel barrier. As tunneling rates are ex- 
ponentially sensitive to tunnel barrier length, we would 
not expect to measure significant displacements from the 
mid-point, since the longer tunnel barrier will limit trans- 
port, and quickly reduce the current below our noise 
floor. This is consistent with the triangulated positions 
that are estimated to have displacements of 5 nm or less 
from the center, for tunnel barriers of 40 — 50 nm total 
length. 

The lateral position is constrained by the relative ca- 
pacitive coupling of the wire vs. the QPC gate. We find 
a much stronger than expected coupling from the QPC 
gate, especially for peaks A and B, where it dominates 
the coupling from the wire gate by a factor of 5 or 10 
times. This ratio pushes the triangulated position far 
from the wire gate, despite its large area. For peaks A 
and B, the simulated positions which best fit the data are 
found below the QPC gate. Even the incomplete set of 
data for peak W allows us to determine its lateral posi- 
tion, though with a larger error bar. Due to the changing 
length of the tunnel barriers as a fuction of lateral posi- 
tion, i.e. the distance between source and drain, easily 



seen in[4J the magnitude of Cs/C and Cd/C also depend 
on lateral position, providing a position for peak W. 

The lateral position of all the resonances are also forced 
closer to or further under the QPC gate in order to pro- 
duce a sufficiently small Cs/C ratio. This is driven 
by the gap between the 2DEGs, which expands out as 
the QPC gate widens. The lateral position below the 
QPC is, therefore, influenced significantly by relatively 
small changes in the edge of the 2DEG, which is one 
of the primary sources of lateral position uncertainty 
under the QPC since the 2DEG position is only esti- 
mated through the semi-classical calculation and the as- 
sumed MIT metallic boundary cut-off. However, even if 
the 2DEG is more depleted than indicated by the semi- 
classical calculation, leading to wider gaps, the donor 
position would still necessarily be under the QPC gate 
because of the Cwire/Cqpc ratio. 

The magnitude of confidence in position provided by 
this triangulation method is sensitive to several factors 
including uncertainty in exact dimensions (e.g., imper- 
fect knowledge about device dimensions) discussed earlier 
as well as three other underlying assumptions discussed 
now. First, the calculation assumes that the Sb elec- 
tron wavefunction radius is 1.8 nm, which implies that 
the electron wavefunction is bulk like. Donors near the 
surface in the presence of vertical electric fields are pre- 
dicted to hybridize to a bound surface state^, which can 
considerably increase the size of the electron distribution 
and shift it to the surface. Second, the calculations as- 
sume the radius of a single electron occupation Sb site, 
although it is possible that the resonance corresponds to 
aA^ = lto7V = 2 (D-) transition. The size of the D- 
wavefunction is considerably larger, probably about 3 
nm radius. Both of these cases represent circumstances 
where the calculations have underestimated the size of 
electron wavefunction. Simulations indicate that an in- 
crease of radial size of order of 2 - 3 does not significantly 
change the position. This highlights an insensitivity of 
using capacitance ratio triangulation to absolute size if 
it is relatively small compared to the tunnel barrier gap 
size. If the total capacitance can be obtained through a 
charging energy measurement, then this also allows size 
determination, which was not available in these cases. 

Triangulation using the measured capacitance ratio 
leads to the conclusion that the only way that resonances 
A or B can be assigned to transport through the D+ to 
DO transition of antimony sites is that the Sb is located 
below the poly gates. It is unlikely that Sb is at these 
locations because the poly is 200 nm further masking 
this region, in addition to the PMMA, from direct im- 
plant and reasonable diffusion lengths for the thermal 
budget are small enough to make it very unlikely that 
Sb is in this region. Defects near the Si/Si02 interface 
are a likely source of a confining potential. Gate ox- 
ides are known to exhibit point like positive charge de- 
fects near the silicon/oxide interface, which would pro- 
duce relatively deep electrostatic confining wells^. Best 
gate oxides are reported to have fixed charge densities on 



the order of 10^^ /cm^, making it probable that over the 
length and width of the QPC gate there will be 1 or 2 de- 
fects. Faster interface state densities, Dit, can be smaller 
densities and are therefore less of a concern. Although 
the confinement of the electron will not be the same as 
that of antimony, 1.8 nm Bohr radius, calculations of 
slightly larger metallic discs, to simulate surface disorder 
dots, are not inconsistent with the capacitance ratios. A 
less likely possible source of the resonance could be ar- 
senic dopants that diffused through the poly and gate ox- 
ide. Simulations using a commercially available process 
simulator (TSUPREM4) indicate that the As concentra- 
tion is well below the boron background concentration of 
r^ 10^^/cm^. If true, it is unlikely that these resonances 
come from arsenic. Implant and diffusion of arsenic is rel- 
atively well understood through years of calibrated mod- 
eling developed for the CMOS industry making simulated 
projections of dopant depth fairly reliable. 



V. SUMMARY 

A tunnel barrier in a silicon split-gate point contact 
was implanted with a small number of Sb donors and 
characterized electrically at T ::^ 300 mK. The device 
fabrication is self-aligned using poly depletion gates. The 
self-aligned process fiow represents a first time proof of 
principle path towards significantly better lateral control 
of donor placement within tunnel barriers and is com- 
patible with single ion detectors^^. The implanted tunnel 
barrier shows considerably more single electron transport 
resonances than an unimplanted, control, case. The ratio 
of the sum of capacitances to the gate capacitance was 
characterized for all available gates for three sharp reso- 
nances. Triangulation of the positions of the resonance is 
made possible through the geometric dependence of the 
capacitances and the positioning of the conducting gates, 
which are configured so that all three spatial directions 
are probed. A semi-classical capacitance model was used 
to bound likely regions that a 1.8 nm radius Sb elec- 
tron would be consistent with the measured ratios. One 
of the resonances is consistent with transport through a 
donor within the implanted constriction. The triangu- 
lation method clearly indicates, on the other hand, that 
the other two resonances, if consistent with a small con- 
fined electron of order of 1.8 nm radius, would have lo- 
cated under the QPC gate. Oxide related defects are 
a likely candidate that might facilitate transport under 
the poly gate although transport through adventitious 
background donors or other impurities cannot be ruled 
out. The triangulation method therefore provides critical 
guidance to distinguish between likely defects and high 
probability candidates related to single electron trans- 
port through intentionally implanted donors. It is fur- 
thermore a demonstration of a method to generally dis- 
tinguish between intentionally engineered single electron 
potentials in contrast to unintentional resonances in more 
disordered systems. 
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